First Measurement of the Muon Neutrino Charged Current 
Quasielastic Double Differential Cross Section 
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A high-statistics sample of charged-current muon neutrino scattering events collected with the 
MiniBooNE experiment is analyzed to extract the first measurement of the double differential cross 
section ( ^^ dcosfl ) charged-current quasielastic (CCQE) scattering on carbon. This result 
features minimal model dependence and provides the most complete information on this process to 
date. With the assumption of CCQE scattering, the absolute cross section as a function of neutrino 
energy {a[Ei,]) and the single differential cross section {-^^) are extracted to facilitate comparison 
with previous measurements. These quantities may be used to characterize an effective axial- vector 
form factor of the nucleon and to improve the modeling of low-energy neutrino interactions on 
nuclear targets. The results are relevant for experiments searching for neutrino oscillations. 
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I. INTRODUCTION 



Neutrino charged-current (CC) scattering without pi- 
ous in the final state is important to measure and char- 
acterize, and is a critical component in the neutrino os- 
cillation program of the MiniBooNE experiment [ll-Q 
at Fermilab. Most of these events are charged-current 
quasielastic scattering (CCQE) of the muon neutrino on 
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a bound nucleon (z/^ + n fi~ + p). A robust model 
of these interactions is required to support future experi- 
ments such as NOvA Q and T2K Q that are also search- 
ing for — > oscillations. Such experiments will use 
interactions to detect the appearance of any i^e 
resulting from oscillations in the large distance between 
production and detection. Additional use will be made of 
CC interactions to normalize the neutrino content at 
production using a near detector and to search for the dis- 
appearance of i'^ via a far detector. These analyses will 
require all available experimental and theoretical insight 
on the CCQE interaction in the sa 1 GeV energy range 
and on nuclear (carbon, oxygen) targets. While many 
unknown quantities are eliminated in these experiments 
by considering ratios of far to near events, the cancella- 
tion is not complete due to differences in neutrino flux 
and backgrounds in the near and far detectors. Thus, 
in order to permit precision oscillation measurements, it 
is important to have an accurate characterization of the 
CCQE differential cross sections over a wide span of neu- 
trino energies. 

Historically, it has proven difficult to accurately define 
the CCQE cross section and precise measurements have 
been unavailable. The experimental execution and data 
interpretation are non-trivial for several reasons. Neu- 
trino beams typically span a wide energy range thereby 
preventing an incoming energy constraint on the reac- 
tion. The neutrino flux itself is often poorly known, 
hampering normalization of reaction rates. Background 
processes are frequently significant and difficult or impos- 
sible to separate from the CCQE signal, for instance, CC 
pion production combined with pion absorption in the 
nucleus. Further complicating the description, the tar- 
get nucleon is not free but bound in a nuclear target and 
correlations between nucleons may be important. There 
are differing detection strategies employed by different 
experiments, for example, some require detection of the 
final state nucleon and some do not. Finally, the nuclear 
target often differs between experiments, thus making 
comparisons less straightforward. 

The current data on CCQE scattering come from a va- 
riety of experiments operating at differing energies and 
with different nuclei Q- Modeling of this data has been 
consistent from experiment to experiment, yet remains 
fairly unsophisticated. Preferred for its simplicity, neu- 
trino CCQE models typically employ a relativistic Fermi 
gas (RFC) model (such as that of Ref. 0) that com- 
bines the bare nucleon physics with a model to account 
for the nucleon binding within the specific nucleus. The 
structure of the nucleon is parametrized with the three 
dominant form factors: two vector, Fi,2(Q^), and one 
axial- vector, Fa{Q'^)- The vector form factors, including 
the (squared four-momentum transfer) dependence, 
are well-determined from electron scattering. The ax- 
ial vector form factor at = is known from neutron 
beta-decay. Neutrino-based CCQE measurements may 
then be interpreted as a measurement of the axial- vector 
mass. Ma, which controls the dependence of Fa, and 



ultimately, the normalization of the predicted cross sec- 
tion. 

This simple, underlying model has led to the situa- 
tion where neutrino CCQE measurements typically pro- 
duce a measurement of Ma independent of neutrino en- 
ergy and target nucleus. The resulting world-average is 
Ma = 1.03 ± 0.02 GeV \^ (a recent summary of the var- 
ious Ma values is provided in Ref. fi(|). It should be 
noted that the data contributing to this world-average 
are dominated by higher precision bubble chamber ex- 
periments using deuterium as a target. In addition, most 
(but not all) of the Ma values have come from the ob- 
served distribution of CCQE events in rather than 
from an overall normalization of the event yield. 

Several experiments have recently reported new results 
on CCQE scattering from high-statistics data samples 
with intense, well-understood neutrino beams. The NO- 
MAD experiment extracted a CCQE cross section and 
Ma from data taken on carbon in the energy range, 
3 < < 100 GeV They observe an Ma value 

and cross section consistent with the prior world-average. 
However, data at lower neutrino energies using carbon or 
oxygen as a target have indicated, through Q^-shape fits, 
a somewhat larger value for Ma (by « 10-30%) [ill-[ll]. 
The SciBooNE experiment has recently reported a pre- 
liminary extraction of the total CCQE cross section on 
carbon that is consistent with these results [3]. To add 
to this, the MINOS experiment has reported a prelimi- 
nary result on Ma extracted from CCQE data collected 
on iron in a similar energy range as NOMAD |15j]. This 
result for Ma is also larger than the world-average. None 
of these experiments has reported differential cross sec- 
tions. 

The CCQE cross section predictions resulting from 
these differing measured values for Ma are too unreli- 
able for use by future oscillation experiments, and the 
collection of existing results remains puzzling. Perhaps 
the currently employed model of the CCQE process is 
too naive and the physics of the bound nucleons can al- 
ter the extracted Ma- There is currently large theoreti- 
cal interest in this area [l^, [l3| and a solution may well 
emerge. Regardless, if the experimental results hold up, 
it is clear that improved measurements will be required 
to sort out the situation. These measurements will need 
to go beyond simply extracting an Ma value as the data 
are already indicating that this strategy is too simplis- 
tic. Experiments should advance to providing cross sec- 
tions, differential where possible, that are correctly nor- 
malized with a predicted neutrino flux (not normalized 
to a different reaction channel in the same data). Experi- 
ments should also strive to reduce the model-dependence 
of their reported results. To this end, we present differ- 
ential cross sections in muon kinematics extracted from 
CCQE events collected in MiniBooNE. 

MiniBooNE has accumulated the world's largest sam- 
ple of I/;, CCQE events (« 150, 000) in the 1 GeV re- 
gion. Using this high-statistics and low-background event 
sample, we report the first measurement of an absolute 
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FIG. 1: (color online) Schematic overview of the MiniBooNE 
experiment including the Booster Neutrino Beamline and 
MiniBooNE detector. 



Vfj, CCQE double differential cross section, the main re- 
sult of this work. In addition, CCQE cross sections in 
several other conventional forms are provided. The lay- 
out of the remainder of this paper is as follows. In Sec- 
tion |lll we provide a summary of the MiniBooNE exper- 
iment, including the Booster Neutrino Beamline (BNB) 
and the MiniBooNE detector. We detail the neutrino 
interaction model used to describe the signal and back- 
ground in Section Hill The CCQE selection and analysis 
strategy is outlined in Section HVl Finally, in Section IVl 
we report the MiniBooNE flux-intcgratcd CCQE double 
differential cross section ( -pf-^-^— 3- ) , the flux-integrated 

^ al f^a cos c*^ ' ' ^ 

CCQE single differential cross section ( J.^ ), and the 

flux-unfolded CCQE cross section as a function of en- 
ergy {a[E^^-^^'~^]). To facilitate comparison with up- 
dated model predictions [l6|,ll7|, we provide the predicted 
MiniBooNE neutrino fluxes and measured cross section 
values in tabular form in the Appendix. 



II. MINIBOONE EXPERIMENT 

A. Neutrino beamline and flux 

The Booster Neutrino Beamline (BNB) consists of 
three major components as shown in Figure [TJ a primary 
proton beam, a secondary meson beam, and a tertiary 
neutrino beam. Protons are accelerated to 8 GeV kinetic 
energy in the Fermilab Booster synchrotron and then 
fast-extracted in 1.6 /is "spills" to the BNB. These pri- 
mary protons impinge on a 1.75 interaction- length beryl- 
lium target centered in a magnetic focusing horn. The 
secondary mesons that are produced are then focused by 
a toroidal magnetic field which serves to direct the result- 
ing beam of tertiary neutrinos towards the downstream 
detector. The neutrino flux is calculated at the detector 
with a GEANT4-based IS] simulation which takes into 
account proton transport to the target, interactions of 
protons in the target, production of mesons in the p-Be 
process, and transport of resulting particles through the 
horn and decay volume. A full description of the calcula- 
tion with associated uncertainties is provided in Ref. [l9| . 
MiniBooNE neutrino data is not used in any way to ob- 
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FIG. 2: (color online) Predicted flux at the MiniBooNE 
detector (a) along with the fractional uncertainties grouped 
into various contributions (b). The integrated flux is 5.16 x 
10"'° //^/POT/cm^ (0 < < 3 GeV) with a mean energy 
of 788 MeV. Numerical values corresponding to the top plot 
are provided in Table FVl in the Appendix. 



tain the flux prediction. The resulting i/^ flux is shown as 
a function of neutrino energy in Figure [2] along with its 
predicted uncertainty. These values are tabulated in Ta- 
ble|V]in the Appendix. The flux has an average energy 
(over 0<E^ <3 GeV) of 788 MeV and comprises 93.6% 
of the total flux of neutrinos at the MiniBooNE detector. 
There is a 5.9% (0.5%) contamination of Df^ {veTi^e)] all 
events from these {non-v^) neutrino types are treated as 
background in this measurement (Section lIVPp . 

The largest error on the predicted neutrino flux re- 
sults from the uncertainty of pion production in the ini- 
tial p— Be process in the target as the simulation pre- 
dicts that 96.7% of muon neutrinos in the BNB are pro- 
duced via TT^" decay. The meson production model in the 
neutrino beam simulation fl9( relies on external hadron 
production measurements. Those of the HARP experi- 
ment [l^l are the most relevant as they measure the tt^ 
differential cross section in p— Be interactions at the same 
proton energy and on the same target material as Mini- 
BooNE. The uncertainty in 7r+ production is determined 
from spline fits to the HARP 7r+ double differential cross 
section data [l^. The spline-fit procedure more accu- 
rately quantifies the uncertainty in the underlying data, 
removing unnecessary sources of error resulting from an 
inadequate parameterization [2l[ of the HARP data. The 
HARP data used was that from a thin (5% interaction 
length) beryllium target run [l^l- While that data pro- 
vides a valuable constraint on the BNB flux prediction, 
additional uncertainties resulting from thick target ef- 
fects (secondary re-scattering of protons and pions) are 
included through the BNB flux simulation. 

The resulting tt"*" production uncertainty is ~ 5% at 
the peak of the flux distribution increasing significantly 
at high and low neutrino energies. There is a small con- 
tribution to the fiux error from the uncertainty in kaon 
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production which is significant only for > 2.0 GeV. 
Other major contributions to the flux error include un- 
certainties on the number of protons on target (POT), 
hadron interactions in the target, and the horn mag- 
netic field. These are grouped as the "beam" component 
shown together with the aforementioned pion and kaon 
production uncertainties in Figure [^b. All flux errors are 
modeled through variations in the simulation and result 
in a total error of w 9% at the peak of the flux. In 
practice, a complete error matrix is calculated in bins of 
neutrino energy that includes correlations between bins. 
This matrix is used to propagate the flux uncertainties 
to the final quantities used to extract the cross section 
results reported here. 



B. MiniBooNE detector 

The MiniBooNE detector (shown schematically in Fig- 
ure [T]) is located 541 m downstream of the neutrino pro- 
duction target and consists of a spherical steel tank of 
610 cm inner radius filled with 818 tons of Marcol? light 
mineral oil (CH2) with a density of 0.845 g/cm^. The 
volume of the tank is separated into an inner and an 
outer region via an optical barrier located at a radius 
of 574.6 cm. The inner and outer regions are only sep- 
arated optically, the oil is the same in each. The in- 
dex of refraction of the oil is 1.47, yielding a Cherenkov 
threshold for particles with /3 > 0.68. The mineral oil is 
undoped, that is, no additional scintillation solutes were 
added. However, because of intrinsic impurities in the 
oil, it produces a small amount of scintillation light in 
addition to Cherenkov light in response to energy loss by 
charged particles. 

The inner (signal) region is viewed by an array of 
1280 inward-facing 8- inch photo multiplier tubes (PMTs) 
mounted on the inside of the optical barrier and providing 
11.3% photocathode coverage of the surface of the inner 
tank region. The outer (veto) region is monitored by 240 
pair-mounted PMTs which record the light produced by 
charged particles entering or exiting the detector volume. 

The PMT signals, in response to the light produced 
from charged particles, are routed to custom-built elec- 
tronics modules where they are amplified, discriminated, 
and digitized. These ("QT") modules extract the start 
time and integrated charge from each PMT pulse that is 
above a discriminator threshold of «0.1 photoelectron. 
This unit of data is called a PMT "hit" . The data is 
stored in a temporary buffer until a trigger decision is 
made. The trigger system uses information from the Fer- 
milab accelerator clock signals and PMT multiplicities to 
form physics and calibration trigger signals. The physics 
trigger for this analysis requires only that beam be sent 
to the BNB neutrino production target. When this con- 
dition is satisfied, all PMT-hit data in a 19.2 fj,s window 
starting 5 fj,s before the 1.6 fis beam spill is extracted 
from the QT modules and added to the data stream. This 
readout strategy collects all PMT data (with no multi- 



plicity threshold) from beam-induced neutrino events as 
well as any muon decays that occur with a characteristic 
time of 2 ^s after the neutrino interaction. 

The data within the 19.2 fj.s readout window is ex- 
amined at the analysis stage to organize the hits into 
temporal clusters or "subevents" . A subevent is any 
group of at least 10 hits (from inner or outer PMTs) 
where no two consecutive hits are separated in time by 
more than 10 ns. These subevents may then be analyzed 
separately to extract further information such as energy 
and position. With this scheme, muon-decay electrons 
or positrons may be identified and separated from the 
primary neutrino event. 

The MiniBooNE detector is calibrated via the light 
from a pulsed laser source, cosmic muons, and muon- 
decay electrons. The laser calibration system consists of 
a pulsed diode laser injecting light via optical fibers into 
four 10 cm dispersion flasks located at various depths in 
the main detector volume. In addition, one bare fiber is 
installed at the top of the tank. This system is used to 
quantify the charge and time response of the PMTs and 
allows for an in situ measurement of the oil attenuation 
length and light scattering properties. 

The muon calibration system consists of a two-layer 
scintillation hodoscope located above the detector tank 
combined with seven 5-cm cubic scintillators deployed 
at various locations within the tank (near the vertical 
axis) and used to tag stopped muons. With this system, 
stopping cosmic ray muons of energies ranging from 20 
to 800 MeV are tracked through the detector enabling an 
energy calibration via the known range of the stopping 
muons. Using this data, the energy (angular) resolution 
for reconstructed muons in MiniBooNE is measured to 
be 12% (5.4°) at 100 MeV, improving to 3.4% (1.0°) at 
800 MeV. 

A separate large sample of stopping cosmic muons is 
obtained via a dedicated calibration trigger that requires 
the signature of an incoming muon and its decay electron. 
This sample allows for the calibration and measurement 
of the detector response to muon-decay electrons. 

More details about the oil medium, detector struc- 
ture, PMTs, electronics, and calibration are available in 
Ref. di. 



C. Detector simulation 

A GEANT3-based |23| program is used to simulate the 
response of the detector to neutrino interactions. This 
simulation is used to determine the detection efficiency 
for CCQE events, the probability for accepting back- 
ground events, and the error on relevant observables due 
to uncertainties in the detector response. 

The entire geometry of the MiniBooNE detector is 
modeled including the detector tank and all inner com- 
ponents. The major components of the detector hous- 
ing structure are modeled as is the surrounding environ- 
ment. Standard GEANt3 particle propagation and decay 
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routines are utilized with some changes made to better 
simulate tt^ decay, ^ decay, and /z~ capture on carbon. 
The latter process is especially important for the back- 
ground estimation in the analysis reported here. The 
GCALOR ^] hadronic interaction package is used in place 
of the default GFLUKA package in GEANtS because it bet- 
ter reproduces known data on 7r+ absorption and charge 
exchange in the relevant 7r+ energy range (100-500 MeV). 

Some modifications were made to the standard 
GCALOR code to better simulate tt^ processes. Pion 
radiative capture/decay and photonuclear processes are 
of concern for the neutrino oscillation search Q but have 
negligible effect for this analysis. Elastic scattering of 
TT* on carbon is important and was not simulated in the 
standard GCALOR code. A model, guided by the avail- 
able data |25i, [22] , was added to the simulation and yields 
only a small change to the calculated background from 
pion production processes. 

Uncertainties of 35% on 7r+ absorption and 50% on 
TT"*" charge exchange are assigned based on the difference 
between the external data ^5j and the GCALOR pre- 
diction. Note that these errors are relevant for 7r+ propa- 
gation in the detector medium not intranuclear processes 
which are assigned separate uncertainties fSection llll Pp . 

Charged particles propagating through the detector oil 
produce optical photons via Cherenkov radiation and 
scintillation. Optical wavelengths of 250-650 nm are 
treated. The Cherenkov process is modeled with stan- 
dard GEANtS. The scintillation process is modeled with a 
MiniBooNE-specific simulation that creates optical pho- 
tons at a rate proportional to Birk's law-corrected energy 
loss with an emission spectrum determined from dedi- 
cated florescence measurements. Optical photons result- 
ing from these production processes are tracked through 
the detector oil with consideration of scattering, fluores- 
cence, absorption, and reflection from detector surfaces 
(including PMTs). Photons that intersect the PMT sur- 
face (and do not reflect) are modeled with a wavelength 
and incident angle- dependent efficiency. The photon sig- 
nal in each PMT is used together with the known re- 
sponse of the PMT and readout electronics to generate 
simulated data that is then input to the data analysis 
programs. 

The models, associated parameters, and errors im- 
plemented in the detector simulation are determined 
through external measurements of the properties of ma- 
terials as well as internal measurements using data col- 
lected in the MiniBooNE detector. It is particularly 
important to correctly model the optical photon trans- 
port since a typical optical photon travels several me- 
ters before detection. This "optical model" is tuned 
starting from various external (to the detector) measure- 
ments of MiniBooNE mineral oil optical properties, such 
as the refractive index, attenuation length, and fluores- 
cence/scintillation strength. These quantities allow for 
the implementation of various models to describe the op- 
tical photon propagation. The details of the models are 
then further adjusted based on MiniBooNE internal data 



sets such as muon-decay electrons, cosmic muons, and 
laser pulses. In total, 35 optical model parameters are 
adjusted to obtain a good description of the various data 
sets. Values for the uncertainties of these parameters, in- 
cluding correlations among them, are also extracted from 
the data and the effect on the reported observables de- 
termined by running the simulation with adjusted values 
(Section HYH. 

Additional details about the MiniBooNE detector sim- 
ulation and supporting measurements of oil properties 
are available in Refs. 122,. ,27J. 



III. NEUTRINO INTERACTION MODEL 

The MiniBooNE experiment employs the nuance v3 
event generator [28j to estimate neutrino interaction rates 
in the CH2 target medium. The nuance generator con- 
siders all interaction processes possible in the neutrino 
energy region relevant for MiniBooNE. It also enables 
the various processes to be tuned to match the data via 
input parameters or source code changes where neces- 
sary. The NUANCE generator includes the following com- 
ponents: (1) a relativistic Fermi gas (RFC) model for 
CCQE (and NC elastic) scattering on carbon Q; (2) a 
baryonic resonance model for CC/NC single and multi- 
pion production [2S|j_(3) a coherent CC/NC single pion 
production model (4) a deep inelastic scattering 

(DIS) model ^3l|, l32j |: and (5) a final-state interaction 
model to simulate re-interaction of final state hadrons in 
the nuclear medium [2^. Neutrino interactions on both 
free (protons) and bound nucleons (in carbon) are con- 
sidered to model the CII2 detector medium. Further de- 
tails on these models, parameters, and uncertainties are 
provided in the following subsections. 

This event generator is used, after the adjustment of 
parameters to adequately describe the MiniBooNE data, 
to calculate the background contribution to the CCQE 
signal. The CCQE model is used only in the extraction of 
the model parameters in a shape-fit to the Qqe distribu- 
tion (Section |TVC|. The CCQE differential cross section 
measurements do not depend on this model (excepting 
some small dependence due to detector resolution cor- 
rections and CCQE backgrounds). In addition, since 
it is such a large background to CCQE, the CCl7r+ back- 
ground is constrained (outside of the nuance model) to 
reproduce MiniBooNE data (Section lIVCp . A summary 
of interaction channels considered and the NUANCE pre- 
dictions for the neutrino interaction fractions in the Mini- 
BooNE detector are provided in Table [H The final values 
for the predicted backgrounds (after event selection cuts) 
are provided in Table WUl 



A. Charged current quasielastic scattering 

CCQE scattering is the dominant neutrino interaction 
process in MiniBooNE and the subject of this analysis. 
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neutrino process abbreviation reaction fraction (%) 

CC quasielastic CCQE i/-^ + n — >■ /j^ + p 39 

NC elastic NCE i/^ +p{n) + p{n) 16 

CC 171+ production CCl7r+ + p{n) ^ /i" + 7r+ + p{n) 25 

CC Itt" production CCItt" v^, + n ^ + n° + p 4 

NC l7r='= production NClTr^l^ +p{n) u^, + 7r+(7r") + n(p) 4 

NC Ivr" production NCItt" v,, +p{n) ^ + tt" + p{n) 8 

multi pion production, DIS, etc. other +p(n) — >■ p^ + A^tt* + X, etc. 4 



TABLE I: Event type nomenclature and NUANCE-predicted i/f^ event fractions for MiniBooNE integrated over the predicted 
flux in neutrino mode before selection cuts. For the pion production channels, indirect production (through resonance states) 
and direct production (through coherent processes) are included. (CC=charged-current, NC=neutral-current). 



This process is defined as the charge-changing scattering 
of a neutrino from a single nucleon with no other parti- 
cles produced and it is simulated with the RFG model [1] 
with several modifications. A dipole form is used for the 
axial form factor with an adjustable axial mass, Ma- An 
empirical Pauli-blocking parameter, k, is introduced 11 1 
to allow for an extra degree of freedom that is impor- 
tant to describe the MiniBooNE data at low momentum 
transfer. This parameter is a simple scaling of the lower 
bound, £'io, of the nucleon energy integral of Ref. [8] via 



E\o = k{^Pp + M^-uj + Eb), where pf, Mp, uj, and Eb 
are the Fermi momentum, nucleon mass, energy transfer, 
and binding energy, respectively. When k > 1, the Pauli- 
blocking of final-state nucleons is increased which reduces 
the cross section at low momentum transfer. 

A parametrization (33l | is used to describe the non- 
dipole behavior of the Dirac and Pauli form factors. Al- 
though the contribution is small, the pseudoscalar form 
factor, derived from partial conservation of the axial vec- 
tor current (PCAC) is also included [s^]. The scalar and 
axial tensor form factors are set to zero as implied from 
G-parity conservation. The Fermi momentum and bind- 
ing energy for carbon are set to 220 ± 30 MeV/c and 
34 ±9 MeV, resp ectively, as extracted from electron scat- 
tering data [3a] taking account of the purely isovector 
character of CCQE. 

The parameters Ma and k were initially extracted 
from MiniBooNE CCQE data in a prior analysis [ll| . 
and were determined to be Mf = 1.23 ± 0.20 GeV and 
K=1.019±0.011. While not the main result of this pa- 
per, this exercise is repeated after explicitly measuring 
the CCItt"'" background from MiniBooNE data and is 
described in Sect IIV CI The superscript "eff " on AIa was 
introduced to allow for the possibility that the axial mass 
measured from scattering on nucleons bound in carbon 
may be different from the "bare-nucleon" axial mass that 
appears within the neutrino model. The use of this nota- 
tion is continued in this work. The uncertainties in these 
CCQE model parameters do not propagate to the errors 
on the measured cross sections for this channel. Neu- 
tral current elastic (NCE) scattering is described with 
the same model as that for the CCQE interaction with 
the replacement of appropriate form factors [ssj to de- 



scribe the NC coupling to the nucleon. The uncertainty 
from the NCE model parameters on the CCQE results is 
negligible due to the small background contribution from 
this channel (Table Hill . 



B. Resonance interactions 

The primary source of single pion production for Mini- 
BooNE is predicted to be baryonic resonance production 
and decay, such as. 



TT +71 , TT + p. 



The NUANCE model employs the relativistic harmonic os- 
cillator quark model of Ref. [l^, [s^ • The pion angular 
distribution due to the spin structure of the resonance 
states is additionally taken into account. In total, 18 
resonances and their interferences are simulated in re- 
actions with invariant mass W < 2 GeV, however, the 
A(1232) resonance dominates at this energy scale. For 
reactions on bound nucleons, an RFG model is employed 
with a uniform Fermi momentum and constant binding 
energy. In-medium effects on the width of resonances 
are not considered explicitly, however, final-state inter- 
actions can produce an effective change in these widths. 
An axial mass of M\'' = 1.10 ± 0.27 GeV, set by tuning 
to available data, is used for this channel. 

Multi-pion processes are considered in the nuance 
simulation with M^'^ = 1.30 ± 0.52 GeV. This parame- 
ter was set (together with M^) so that the simulation 
reproduces inclusive CC data. The contribution of this 
channel to the CCQE background is small and the un- 
certainty is negligible in the final errors. 

The CCl7r+ channel is the largest background contri- 
bution to the CCQE signal and the uncertainty from the 
model prediction alone is substantial [s^. However, the 
experimental signature of the CCl7r+ reaction in Mini- 
BooNE is distinct and the efficiency is large. So, in order 
to reduce uncertainty stemming from the CCl7r+ model 
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prediction, a measurement of the CCl7r+ background is 
performed as part of the CCQE analysis (Section IIV Cp . 



C. Coherent pion production 

Pions are also produced in the CC and NC coherent 
interaction of neutrinos with carbon nuclei, 



A 
A 



^0 



f A 

A . 



In NUANCE, this coherent pion production process is de- 
scribed using the model of Ref. [s^ assuming M'j^^ = 
1.03 ±0.28 GeV [H. 

Coherent scattering is predicted to have distinct fea- 
tures in the angular distributions of both the final- 
state muons and pions. Both the K2K [39] and Mini- 
BooNE [i^ experiments have measured the fraction of 
pions produced coherently in 1 GeV neutrino interac- 
tions. K2K measured a rate for coherent CCl7r+ produc- 
tion consistent with zero and set an upper limit. Mini- 
BooNE measured a non-zero value for coherent NCItt*^ 
production albeit « 35% smaller than the model pre- 
diction m, [S^. The latest resuh from the SciBooNE 
experiment is consistent with the K2K measurement for 
CCl7r+ coherent production [4l|. 

Because of the current discrepancy between CC and 
NC coherent pion measurements and the variation in 
model predictions at low energy, the prediction is reduced 
to 50% of the default value in nuance [H, H^I and as- 
signed a 100% uncertainty. This choice spans the current 
results from relevant experiments and existing theoretical 
predictions. 



D. Final State Interactions 

In the NUANCE simulation, neutrino interactions on 
nucleons are modeled using the impulse approximation 
which assumes the interaction occurs instantaneously on 
independent nucleons. The binding of nucleons within 
carbon is treated within the RFC model, however, any 
nucleon-nucleon correlation effects are not. The final 
state hadrons may interact within the nucleus as they 
exit. They are propagated through the ^^C nucleus with 
a known, radially-dependent nucleon density distribu- 
tion 42] and may undergo final-state interactions (FSI). 
These are simulated by calculating interaction probabil- 
ities for the possible processes in 0.3 fm steps until the 
particles leave the ss 2.5 fm-radius spherical carbon atom. 
The interaction probabilities are derived from external 
TT — N , N — N cross section and angular distribution 
data Q, as well as the nuclear density of carbon. 

To model A absorption in the nucleus, (A -1- ^ 
iV -f N), a constant, energy-independent probability for 
an intranuclear interaction of 20% (10%) is assumed for 
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FIG. 3: A comparison of relevant data [2^ with the nu- 
ance model (solid lines) for intranuclear pion (a) absorption 
and (b) charge exchange as a function of pion kinetic energy. 
The dotted lines show the 25% (30%) systematic error bands 
assumed for the pion absorption (charge-exchange) cross sec- 
tions. 



A+-|-A^, ISP+N (A++-fiV, I\~+N) processes. These val- 
ues were chosen based on comparisons to K2K data [i^l 
and are assigned a 100% uncertainty. After an interac- 
tion, the density distribution and step size are modified 
to prevent an overestimate of these FSI effects [i^ . 

Of these hadron FSI processes, intranuclear pion ab- 
sorption and pion charge exchange, (tt^ ~\- X ^ X' ,tt^ + 
X ^ -K^ + X') are the most important contributions 
to the uncertainty in the CCQE analysis. Pion ab- 
sorption and charge exchange in the detector medium 
are addressed separately in the detector simulation (Sec- 
tion |TrC|). 

A CCl7r+ interaction followed by intranuclear pion ab- 
sorption is effectively indistinguishable from the CCQE 
process in MiniBooNE (they are "CCQE-like", Sec- 
tion II VP I) . An event with intranuclear pion charge ex- 
change is distinguishable in the detector, albeit not with 
100% efficiency. These effects, combined with the high 
rate of CCl7r+ events, results in a significant background 
to the CCQE measurement that must be treated care- 
fully. The model for these intranuclear pion processes 
has been tuned to match the available data 25] in the 
relevant pion energy range. A comparison of the adjusted 
NUANCE model and relevant data is shown in Figure [S] 
A 25% (30%) systematic error in the overall interaction 
cross section is used for the pion absorption (charge- 
exchange) process. 
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FIG. 4: (color online) Schematic illustration of a CCQE inter- 
action in the MiniBooNE detector. The primary Cherenkov 
light from the muon (Cherenkov 1, first subevent) and sub- 
sequent Cherenkov light from the decay-electron (Cherenkov 
2, second subevent) are used to tag the CCQE event. No 
requirements are made on the outgoing proton. 



IV. CCQE MEASUREMENT 

The goal of this measurement is to determine the dou- 
ble differential cross section for the CCQE process on 
carbon, + n ^ fi^ + p, vi^herc the target neutron is 
bound in ^^C. 

The identification of CCQE interactions in the Mini- 
BooNE detector relies solely on the detection of the 
Cherenkov light from the primary (prompt) muon and 
the associated decay-electron. An illustration of this 
process is shown in Figure S) Scintillation light is pro- 
duced by the charged lepton and the recoil proton (or 
nuclear fragments). However, with the reconstruction 
employed here, this light is not separable from the dom- 
inant Cherenkov light. In addition, the proton is typi- 
cally below Cherenkov threshold. These conditions are 
such that the proton is not separable from the charged 
lepton and so no requirement is placed on the recoil pro- 
ton in this analysis. This is to be contrasted with some 
measurements of CCQE interactions that do require the 
observation of a recoil proton for some part of the event 
sample [13, fl^ - fl^ . An advantage of this insensitivity 
to the proton recoil is that the extracted cross sections 
are less dependent on proton final-state model uncertain- 
ties. However, the disadvantage in not detecting the re- 
coil nucleon is that contributions to scattering from other 
nuclear configurations (such as two-nucleon correlations) 
are inseparable. These contributions are, in the strictest 
sense, not CCQE, but counted as such in our experimen- 
tal definition. 

A requirement of low veto activity for the CCQE sam- 
ple ensures that all particles produced in the event stop 
in the main region of the detector. This allows muons 
to be tagged with high efficiency via their characteristic 
electron-decay with r « 2 /is. 

The CCQE interaction, including the muon decay, pro- 



where each line in this equation identifies a subevent 
(Section III Bp . The primary muon is identified with the 
first subevent and the subsequent decay-electron with the 
second subevent. At BNB neutrino energies, neutrino 
interaction events that contain a primary muon predom- 
inantly result from CCQE scattering as can be seen in 
Table H 

The largest background is from CC single-pion pro- 
duction (CCl7r+). A CCl7r+ interaction in the detector 
consists of (with subevents labeled), 



1: i^f^ + p{n) fi + p{n) + n'^ 



2 : 

3 : 



Note that this interaction results in three subevents: the 
primary interaction and two muon decays (the muon de- 
cays can occur in any order). The 7r+ decays immediately 
and light from the prompt decay products contribute to 
the total light in the primary event. These events may 
be removed from the CCQE sample by requiring exactly 
two subevents. This requirement also reduces the back- 
ground from NC processes to an almost negligible level 
because they do not contain muons and thus have only 
one subevent. This simple strategy results in a fairly pure 
sample of CCQE events. However, a significant number 
of CCl7r+ events have only two subevents because one of 
the decay electrons escapes detection: the fj,~ is captured 



on ^^C in the mineral oil (with 8% probability [47[) or 



the 7r+ is absorbed. Additionally, the study of CCItt"' 
events for this analysis has indicated that the prediction 
for the CCl7r+ channel from the nuance event genera- 
tor is not sufficiently accurate for this measurement [s^. 
For these reasons, the CCl7r+ rate is measured using 
a dedicated event sample. This differs from our previ- 
ous strategy [ll| where the default NUANCE-predicted 
CCl7r+ fraction (with no adjustments) was used, with 
generous errors, in fits to the CCQE sample. 

The resulting procedure for selecting the CCQE sam- 
ple and measuring the CCl7r+ background involves the 
following steps: 

1. selection of a "super-sample" of events with a clean 
muon signature to isolate CC events (predomi- 
nantly CCQE and CCl7r+) via analysis cuts; 

2. application of a subevent cut to separate the super- 
sample into CCQE (2-subevents) and CCItt^ (3- 
subevents) samples; 



3. measurement of the CCl7r+ 
sample; 



rate from the CCItt"' 
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4. adjustment of the CCl7r+ model in the event sim- 
ulation to reproduce the measured rate; and 

5. subtraction of this adjusted CCl7r+ background 
(along with other predicted backgrounds) from the 
CCQE signal to produce a a measurement of the 
CCQE interaction cross section. 

The details of this procedure are provided in the following 
subsections. 

In this analysis, the reconstruction of the CCl7r+ sam- 
ple is for the sole purpose of background estimation. Ded- 
icated measurements of the CCl7r+ and CCItt" channels 
in MiniBooNE have been reported elsewhere [48l - l50| in- 
cluding detailed reconstruction of the tt"*" and n'^ kine- 
matics. 



A. Event reconstruction 

For this analysis, it is crucial to identify and measure 
the muon in the CC interaction. This is accomplished 
with an "extended-track" reconstruction algorithm fEl] 
which uses the charge and time information from all PMT 
hits in the first subevent to form a likelihood that is maxi- 
mized to determine the best single track hypothesis quan- 
tified by the track starting point, starting time, direction, 
and kinetic energy. This is performed with both a muon 
and electron particle hypothesis from which a (log) like- 
lihood ratio is formed to enable particle identification. 

The muon kinetic energy, T^, and muon scattering 
angle, 0^, are extracted from the track reconstruction 
assuming a muon hypothesis. These are used to form 
the fundamental observable reported here, the double- 
differential cross section. For additional reported observ- 
ables, the reconstructed neutrino energy E^^ and recon- 
structed four-momentum transfer Qq^ are obtained via, 

j^QE ^ 2{M[,)E^-{{M'^f+ml-Ml) 

2 • [(M;) -E^ + ^El - ml cos 9^] ' 

Qqe = -K + 2i??''(i?M - ^El-mlcose^), (2) 

where E^ = + is the total muon energy and Af„, 
Mp, are the neutron, proton, and muon masses. The 
adjusted neutron mass, M'^ = M„ — Eb^ depends on the 
binding energy (or more carefully stated, the separation 
energy) in carbon, Eb, which for this analysis is set to 
34 ± 9 MeV. 

The subscript, '^QE^\ on these reconstructed quanti- 
ties is to call attention to these specific definitions and to 
distinguish them from quantities obtained in other ways 
such as fits to the underlying true kinematic quantities. 
These are kinematic definitions that assume the initial 
nucleon (neutron) is at rest and the interaction is CCQE 
( "QE assumption" ) . While these quantities certainly dif- 
fer from the underlying true quantities, they are well- 
defined, unambiguous, and easily reproduced by others. 



B. CCQE and CCl7r+ event selection 

The CCQE and CCl7r+ candidate events are selected 
for this analysis and separated with a sequence of cuts 
summarized in Table HIl 

The first five cuts are designed to efficiently select a 
high-purity sample of CCQE and CCl7r+ events. Cut 1 
rejects events with incoming particles such as cosmic rays 
or neutrino-induced events produced in the surrounding 
material. It also eliminates events where any of the neu- 
trino interaction products escape the main detector vol- 
ume. This is important for an accurate muon energy 
measurement and to avoid missing muon decays which 
leads to higher backgrounds. Cut 1 does reduce the effi- 
ciency substantially (Tab. HI)) , however, it is necessary to 
reduce background (together with the subsequent cuts). 
Cut 2 requires that the primary (muon) is in-time with 
the BNB spill window. Cut 3 ensures that the recon- 
structed primary muon vertex is located within a fidu- 
cial region in the main detector volume sufficiently far 
from the PMTs for accurate reconstruction. Cut 4 pro- 
vides a minimum muon kinetic energy for reliable recon- 
struction and reduces backgrounds from beam-unrelated 
muon-decay electrons. 

Cut 5 requires that the candidate primary muon is 
better fit as a muon than as an electron. Misrecon- 
structed and multi-particle events tend to prefer the elec- 
tron hypothesis so this cut reduces such contamination. 
This also substantially reduces the efficiency for selecting 
CCl7r+ events as can be seen in Figure [S] where the /i/e 
log-likelihood ratio distribution is shown for each of the 
2- and 3-subevent samples. This bias is intended as it 
selects a sample of CCItt"*" with muon kinematics more 
closely matched to those CCl7r+ that are background 
to the CCQE sample. As is shown in Fig. [5l data and 
Monte Carlo simulation (MC) agree fairly well to within 
the detector errors. The log-likelihood ratio distribution 
is quite sensitive to details of an event such as scintil- 
lation from hadron recoil via the PMT charge and time 
information [5l| . The data-MC difference in the number 
of events passing Cut 5 in both the 2- and 3-subevent 
samples is covered by the full systematic errors consid- 
ered in this analysis. 

Cut 6 separates the samples into CCQE (2 subevents) 
and CCl7r+ (3 subevents) candidates. For this analysis, 
the second and third subevents are required to contain at 
least 20 tank hits to reduce the probability of accidental 
coincidences with the initial neutrino interaction (first 
subevent). This requirement reduces the efficiency for 
identifying the muon-decay electron by « 3%. 

Cut 7 utilizes the muon-electron vertex distance, the 
measured separation between the reconstructed muon 
and electron vertices. This cut requires that the decay- 
electron is correctly associated with the primary muon 
and is applied to the CCQE (2-subevent) sample only. 
This eliminates many CCItt^ events where the second 
subevent is a decay-positron from the tt"*" decay chain and 
not the electron from the decay of the primary muon. 
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CCQE 


cut # 


description 


efEc.(%) purity (%) 


1 


all subevents, # of veto hits < 6 


54.8 


36.8 


2 


1st subevent, event time window, 4400 < T(ns) < 6400 


54.3 


36.8 


3 


1st subevent, reconstructed vertex radius < 500 cm 


45.0 


37.4 


4 


1st subevent, kinetic energy > 200 MeV 


39.7 


46.3 


5 


1st subevent, /i/e log-likelihood ratio > 0.0 


36.0 


62.3 


6 


# total subevents = 2 for CCQE (= 3 for CCl7r+) 


29.1 


71.0 


7 


(CCQE-only) 1st subevent, fi — e vertex distance > 100 cm and 








^ - e vertex distance > (500 x Tp(GeV) - 100) cm 


26.6 


77.0 



TABLE II: List of cuts for the CCQE and CCItt"^ event selections. The predicted efficiency and purity values are for the CCQE 
signal normalized to all CCQE events with a reconstructed vertex radius, r < 550 cm. 




-0.1 -0.08 -0.06 -0.04 -0.02 0.02 0.04 0.06 0.08 logL(n/e) 

FIG. 5: (color online) Distributions of the fi/e log-likelihood 
ratio for the (a) 2 and (b) 3 subevent samples. Data and 
Monte Carlo simulation (MC) are shown along with the in- 
dividual MC contributions from CCQE, CCItt^, and other 
channels. The lines with arrows indicate events selected by 
the muon-electron likelihood ratio cut. 



The distributions of the muon-electron vertex distance 
for the major Monte Carlo channels and for data are 
shown in Figure [6l after the application of cuts 1-6 for 
events with 2 subevents. 

As shown in Table HIl the efficiency for finding CCQE 
events with a reconstructed vertex radius, r < 550 cm, is 
26.6%. An r < 550 cm volume is used for normalization 
in the cross section calculations to correctly account for 
events that pass all cuts but have a true vertex with 
r > 500 cm. Normalizing to events with true vertex of 
r < 500 cm yields an efficiency of 35%. 



C. CCItt background measurement 

After the selection of the CCQE and CCl7r+ candi- 
date events (2 and 3 subevent samples, respectively), the 
CCItt"'' background to the CCQE signal is measured by 




T,(GeV) T^(GeV) 



FIG. 6: (color online) Scatter plots of muon-electron vertex 
distance as a function of reconstructed muon kinetic energy 
for MC samples: (a) CCQE, (b) CCl7r+, and (c) other chan- 
nels. The distribution of data is shown in (d). These are 
events with 2-subevents and with cuts 1-6 applied. The lines 
with arrows indicate the region selected by the muon-electron 
vertex distance cut. 



adjusting the weights of the simulated CClTr"*" events to 
achieve data-MC agreement in the Qqe distribution of 
the 3 subevent sample. The same weighting, applied to 
all simulated CCItt"*" events, then provides an estimate 
of the CCItt"*" background to the CCQE signal. Figure [7] 
shows the data and MC Qqe distributions for the two 
samples before the reweighting of CCItt"*" MC events. 
The 3-subevent sample is predicted to be 90% CCItt"*" 
and shows a large data-MC disagreement in both shape 
and normalization. The kinematic distribution of muons 
in CCItt"*" events is similar in both the 2-and 3-subevent 
samples as can be observed in Fig. [71 This occurs because 
the majority of CCItt"*" events that are background in the 
2-subevent sample are due to muon-capture or pion ab- 
sorption and the reconstruction of the primary event is, 
to a good approximation, independent of this. In addi- 
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FIG. 7: (color online) Distribution of events in Qqe for the 
(a) 2 and (b) 3 subevent samples before the application of 
the CCItt^ background correction. Data and MC samples 
are shown along with the individual MC contributions from 
CCQE, CCItt+j and other channels. In (b), the dashed line 
shows the CCItt"^ reweighting function (with the y-axis scale 
on the right) as determined from the background fit proce- 
dure. 



tion, the fi/e log-likelihood ratio cut (Tab. Ulland Fig. [5]) 
is applied for both the 2- and 3-subevent samples, fur- 
ther ensuring that the CCl7r+ events are the same in 
both samples. 

The CCl7r+ reweighting function (Fig. Uh) is a 4**^- 
order polynomial in Qq^ and is determined from the 
ratio of data to MC in this sample. The 2-subevent 
sample shows good shape agreement between data and 
MC. This is because the event model for CCQE was al- 
ready adjusted to match data in a previous analysis [llj | 
that considered only the shape of the Qq^ distribution. 
That analysis did not consider the overall normalization 
of events. 

In practice, this determination of the CCl7r+ reweight- 
ing is done iteratively as there is some CCQE background 
in the 3 subevent sample. An overall normalization factor 
is calculated for the CCQE sample to achieve data-MC 
agreement in the 2 subevent sample after subtraction of 
the CCItt"'' background. This is then applied to deter- 
mine the CCQE background in the 3 subevent sample. 
The background from other channels is determined from 
the simulation and subtracted. This process converges 
after two iterations. 

This method determines a correction to the CCl7r+ 
rate (as a function Qqe) using data from the 3-subevent 
sample rather than relying strictly on simulation. This 
reweighting is then applied to all simulated CCl7r+ 
events, in particular those that are contained in the 2- 
subevent sample and form most of the background for 
the CCQE measurement. The error on within the 
resonant background model is then set to zero and the 
resulting error on the CCl7r+ background to the CCQE 
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FIG. 8: (color online) Distribution of events in Qqe for the 
(a) 2 and (b) 3 subevent samples after the application of the 
data-based CCItt"^ background constraint and the new CCQE 
model parameters M'^ and k as determined from the CCQE 
fit procedure described in the text. 



signal from CCItt^ production is determined by the co- 
herent TT-production errors and the 7r+ absorption un- 
certainty. The statistical errors in this procedure are 
negligible. Most CCItt"*" events that end up in the 2- 
subevent (CCQE) sample are due to intranuclear tt"*" ab- 
sorption. This process is modeled in the event simulation 
as explained in Sec. IIIIDI and is assigned a 25% uncer- 
tainty. The coherent 7r-production process is modeled as 
described in Sec. lIII Cl aud is assigned a 100% uncertainty. 

With the measured CCl7r+ background incorporated, 
a shape-only fit to the 2-subevent (CCQE) sample is per- 
formed to extract values for the CCQE model parame- 
ters, and k. This exercise is required to have a 
consistent description of the MiniBooNE data within the 
simulation after adjustment of the background. This pro- 
cedure has no effect on the CCQE cross section results 
reported here other than very small corrections to the an- 
tineutrino background subtraction which uses these pa- 
rameters. In this fit, all systematic errors and correla- 
tions are considered. The CCQE simulated sample is 
normalized to have the same number of events as data 
which is the same normalization as determined in the 
CCItt"'" background determination. The Qq^ distribu- 
tions of data from the 2 and 3 subevent samples is shown 
together with the MC calculation in Figure M The MC 
calculations include all the adjustments described in this 
section and agreement with data is good in both samples. 

This shape-only fit to the 2-subevent sample yields the 
adjusted CCQE model parameters, and k, 

Mf = 1.35 ±0.17 GeV/c^ ; 
K = 1.007 ±0.012 ; 
xVdof = 47.0/38 . 

Figure [ni shows the Icr contour regions of this fit together 
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FIG. 9: (Color online). The Icr contour plot for the Mf - k 
fit. The filled star shows the best fit point and Icr contour 
extracted from this work. The open star indicates the best fit 
point and Icr contour from the previous work (TJ]. Two error 
ellipses are shown for the previous work, the larger reflects 
the total uncertainty ultimately included in Ref. The 
circle indicates the world-averaee value for Ma ^■ 



with the results from the original MiniBooNE analy- 
sis [m . The new fit yields different results for both Mf 
and K because of the improved CCl7r+ background es- 
timation method used in this analysis. Note that the 
current result is consistent to within Icr with k — 1, un- 
like the previous MiniBooNE result. This fit actually 
provides no lower bound on k, as the la contour is not 
closed for k < 1. The value for k is quite sensitive to 
the CCl7r+ background at the lowest Qg^ and the back- 
ground in that region has decreased in this analysis. The 
increase in the CCl7r+ background at larger Qg^ val- 
ues has resulted in a larger value for the extracted Mf . 
The previous and current parameter contours are con- 
sistent at the Icr level. Neither this nor the prior anal- 
ysis result is consistent with the world-average Ma of 
1.03 ± 0.02 GeV @, as can be seen in Figure El The 
X^/dof assuming Ma = 1.03 GeV, k = 1 is 67.5/40 cor- 
responding to a probability of w 0.5%. 

The reconstructed four-momentum transfer, Qq^, de- 
pends upon the muon energy as can be seen in Eq. |2j The 
reconstructed muon energy calibration has been checked 
by comparing the measured range of muons determined 
from the muon-electron vertex distance (Fig. ^ with the 
expected muon range determined from the energy pro- 
vided by the reconstruction algorithm, which does not 
use this vertex distance. As an example, a comparison 
of the measured and expected muon ranges for muons 
of 400 < Tf, < 500 MeV is shown in Fig. [TO] for both 
data and simulation. The agreement is good for all muon 
energies and verifies the energy calibration to 2%, well 
within the errors calculated by the simulation. This also 
shows that any light produced by hadronic particles in 



FIG. 10: (Color online). Distributions, for both data and 
simulation, of the measured muon range for muons with 400 < 
< 500 MeV. The subtracted offset, (500 x Tfj.{GeV) - 
100) cm, corresponds to the solid diagonal line in Fig. [6] 
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FIG. 11: Ratio of CCQE data to simulation as a function 
of reconstructed muon kinetic energy and angle. The simu- 
lation includes adjustment of the CCItt"'' backgrounds based 
on MiniBooNE data but is that prior to any CCQE model 
tuning (it assumes A'Ia ~ 1.03 GeV, k — 1). The prediction 
has been normalized to the data. If the simulation modeled 
the data perfectly, this ratio would be unity. Contours of con- 
stant El, and are overlaid, and only bins with > 20 events 
in the data are plotted. 



the neutrino interaction (for both CCQE and background 
channels) is adequately simulated and considered in the 
reconstruction. 

A final and more complete check that the simulation 
correctly models the data can be made by examining the 
2-dimensional muon kinetic energy and angle (T^, cos^^) 
distributions. While Figure |8] shows that data is well- 
described in Qqe, the adjusted model may not be ad- 
equate when applied to the (T^,cos6'^) distribution of 
events. This could occur if an adjustment in Qq^ is 
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FIG. 12: Ratio of CCQE data to simulation as a function 
of reconstructed muon kinetic energy and angle. The simula- 
tion includes the adjusted CCItt"'" background prediction and 
the new CCQE model parameters {M^ = 1.35, k = 1.007) 
measured from MiniBooNE data . Compare to Figure [TT] 



hiding an incorrect neutrino energy distribution. The ra- 
tio of data to simulation in (T^, cos 6*^) after correction of 
the CCl7r+ background and before the new CCQE model 
parameters are applied is shown in Figure 1111 The ratio 
after all corrections is shown in Figure [12] and is much 
closer to unity throughout the muon phase space. As can 
be seen in Figure [TTl the regions of constant ratio mainly 
follow lines of constant Qq^ and not E^^ . Also, almost 
no structure remains in Figure [T^ The exception is band 
of w 20% disagreement at E^^ « 0.4 GeV where the er- 
ror on the neutrino flux is of that order (Fig. [2}d.) This 
data-simulation agreement provides additional support 
for our procedure of adjusting only the Qqe behavior of 
the model and not the energy distribution of the neutrino 
flux. 



D. Extraction of the cross sections 

With the CCItt''" interaction background prediction 
determined from experimental data (Section IIV Cp and 
the remaining channels predicted from the interaction 
model (Section the cross section for the CCQE in- 
teraction can be extracted. A total of 146070 events 
pass the CCQE selection (Section IIVB[) resulting from 
5.58 X 10^" protons on target (POT) collected between 
August 2002 and December 2005. The efflciency for 
CCQE events passing these cuts is calculated to be 26.6% 
for CCQE events with true vertices within a 550 cm ra- 
dius from the center of the detector tank. The sample 
is estimated to contain 23.0% background events. These 
numbers, together with a breakdown of predicted back- 
grounds, are summarized in Tabic Hill 

The background is dominated by CCl7r+ interactions 
which are estimated to be 18.4% of the CCQE candidate 



integrated protons on target 




5.58 X 10^° 


energy-integrated v^j, flux 


2.88 X 


10" u^/cm^ 


CCQE candidate events 




146070 


CCQE efliciency {R < 550 cm) 




26.6% 


background channel 


events 


fraction 


NCE 


45 


< 0.1% 


CCl7r+ 


26866 


18.4% 


CCl7r° 


3762 


2.6% 


NCItt* 


535 


0.4% 


NCl7r° 


43 


< 0.1% 


other 


328 


0.2% 


all non-i^^j 


1977 


1.4% 


total background 


33556 


23.0% 



TABLE III; Summary of the final CCQE event sample includ- 
ing a breakdown of the estimated backgrounds from individ- 
ual channels. The fraction is relative to the total measured 
sample. The channel nomenclature is defined in Table H] 



sample. Their predicted distribution in Qq^ is shown 
in Figure EJa) . As can be seen, this background is a 
substantial fraction of the sample in the lowest Qq^ re- 
gion. The majority (52%) of the CCItt"*" background is 
predicted to be events in which the tt"^ is absorbed in 
the initial target nucleus. These are defined as "CCQE- 
like" in that they contain a muon and no pions in the 
final state. The remaining CCItt"*" background consists 
of CC1tt+ events where the tt"*" is absorbed outside of the 
target nucleus (33%), is not identified due to a missed 
muon decay (11%), or undergoes charge exchange in the 
nucleus or detector medium (4%), These last three classes 
of CCItt"*" backgrounds are not considered CCQE-like. 
All CCItt"*" background events, including those that are 
CCQE-like, are subtracted from the data to obtain the 
final CCQE cross section results. However, to facilitate 
examination of the model used for these processes, the 
effective cross section for CCQE-like background events 
is separately reported in the Appendix. 

In this analysis, the small contamination of P^, z/e, and 
Df, interactions are treated as background and are sub- 
tracted from the data based on their MC prediction (see 
Tab. imp . The majority of these are CCQE inter- 
actions. The same and k as measured in the 
CCQE sample are used to predict non-v^ CCQE events. 
These parameters have been shown to adequately repro- 
duce the MiriiBooNE CCQE data collected in antineu- 
trino mode [52]. 

To extract differential cross sections in muon kine- 
matic variables, the reconstructed kinematics are cor- 
rected for detector-specific effects. A correction proce- 
dure is implemented using an "unfolding" process based 
on the detector simulation. We employ an "iterative 
Bayesian" method [s^l to avoid the problem of amplifi- 
cation of statistical fluctuations common in the "inverse 
response matrix" method [sj]. A disadvantage to the it- 
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erative Bayesian method is that the result depends on 
the initial CCQE model assumptions (the "prior" prob- 
ability). However, this problem is addressed by an iter- 
ative method that uses the extracted signal distribution 
to correct the predicted distributions and repeating this 
procedure. In practice, the simulation was already tuned 
to reproduce the data based on previous work |ll| . and 
the result from the first iteration shows satisfactory con- 
vergence. The systematic uncertainty in this procedure is 
determined from the difference between the initial (0th) 
and final (1st) iterations of the algorithm and by examin- 
ing the dependence of the final values on the initial model 
assumptions. 

The various correction and normalization factors can 
then be brought together in a single expression used to 
extract the CCQE cross section for the ith bin of a par- 
ticular kinematic variable, 



(3) 



where the index j labels the reconstructed bin and i la- 
bels the unfolded (estimate of the true) bin. In this equa- 
tion, Uij is the unfolding matrix, is the efficiency, T is 
the number of neutrons in the fiducial volume, and $ is 
the neutrino flux. This expression is used to obtain the 

,2 

double and single difi^erential cross sections, ,^ ^ „ , 

^ ' dT^d cos 0^ ' 

and -tM—, respectively, after the multiplication of the 

appropriate bin width factors. Note that the choice of 
normalization yields cross sections "per neutron" . Here, 
the flux, $, is a sing le number (2.88 x 10" v^/cxn^) 
and is determined by integrating the BNB flux over 
< -E^ < 3 GeV. Therefore, these differential cross sec- 
tions are "flux-integrated" . 

An additional quantity, the flux-unfolded CCQE cross 
section as a function of neutrino energy, a[E^^'^^'~^\^ is 
extracted from this same expression with the replacement 
of the total flux, with the flux in a particular neutrino 
energy bin, $i. The unfolding procedure is used to cor- 
rect the data from bins of reconstructed neutrino energy, 
E^^ , (using Eq. [T|) to an estimate of the true neutrino 
energy, E^^-^^'^ . It is important to note that, unlike 
the differential cross sections. 



d'a 



and 



the 



calculation of this cross section relies on the interaction 
model to connect E^^ to E^^^^^^. The superscript 
''''RFG^' indicates the interaction model assumed in the 
unfolding process ^ . This procedure introduces a model 
dependence into this cross section, however, this method 
is consistent with that commonly used by experiments 
reporting a CCQE cross section as a function of neu- 
trino energy. This model-dependence should be consid- 
ered when comparing measurements of this quantity from 
different experiments. 



E. Error analysis 

The errors on the measured cross sections result from 
uncertainties in the neutrino flux, background estimates, 
detector response, and unfolding procedure. To propa- 
gate these error sources, a "multisim" method [s^ is used 
to calculate the errors on the final quantities by varying 
parameters in separate simulations. This method pro- 
duces an error matrix, Vij , for each reported distribution 
that can then be used to extract the error on each bin 
(fi — \fVii) and the correlations between quantities of 
different bins. 

The error matrix is calculated by generating a large 
number of simulated data sets with different parame- 
ter excursions, based on the estimated Icr uncertainties 
in those parameters and the correlations between them. 
The error matrix for a particular distribution is then cal- 
culated from these M data sets. 



1 



M 



s=l 



(4) 



Here, Q\ is the quantity of interest in the i*'' bin from the 
s*'* simulation data set and Qi is the "best" estimate of 
the parameters. The quantities of interest Qi could be, 
for instance, the number of events in each bin of Qq^ or 
the calculated cross section in each kinematic bin (Eq.[3]). 

In practice, the errors are classified into four major 
contributions: the neutrino flux, background cross sec- 
tions, detector modeling, and unfolding procedure. The 
parameters within each of these groups are varied inde- 
pendently so the resultant error matrices from the indi- 
vidual groups can be added to form the total error ma- 
trix. For the neutrino flux and background cross section 
uncertainties, a reweighting method is employed which 
removes the difficulty of requiring hundreds of simula- 
tions with adequate statistics. In this method, each neu- 
trino interaction event is given a new weight calculated 
with a particular parameter excursion. This is performed 
considering correlations between parameters and allows 
each generated event to be re-used many times saving 
significant CPU time. The nature of the detector uncer- 
tainties does not allow for this method of error evaluation 
as parameter uncertainties can only be applied as each 
particle or optical photon propagates through the detec- 
tor. Approximately 100 different simulated data sets are 
generated with the detector parameters varied according 
to the estimated Icr errors including correlations. Eq.Uis 
then used to calculate the detector error matrix. The er- 
ror on the unfolding procedure is calculated from the dif- 
ference in final results when using different input model 
assumptions fSection llVPp . The statistical error on data 
is not added explicitly but is included via the statistical 
fluctuations of the simulated data sets (which have the 
same number of events as the data). 

The flnal uncertainties are reported in the following 
sections. The breakdown among the various contribu- 
tions are summarized and discussed in Section [VDI For 
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MiniBooNE data (5Ni.=10.7%) 
MiniBooNE data with shape error 
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FIG. 13: (Color online). Flux- integrated double differential 
cross section per target neutron for the v,j, CCQE process. 
The dark bars indicate the measured values and the surround- 
ing lighter bands show the shape error. The overall normal- 
ization (scale) error is 10.7%. Numerical values are provided 
in Table IVTl in the Appendix. 



simplicity, the full error matrices are not reported for all 
distributions. Instead, the errors are separated into a to- 
tal normalization error, which is an error on the overall 
scale of the cross section, and a "shape error" which con- 
tains the uncertainty that does not factor out into a scale 
error. This allows for a distribution of data to be used 
(e.g. in a model fit) with an overall scale error for un- 
certainties that are completely correlated between bins, 
together with the remaining bin-dependent shape error. 



MiniBooNE data with sliape error 
RFG model (Mf =1.03 GeV,K=1.000) 
RFG model (Mf =1.35 GeV,K=1.007) 
RFG model (Mf =1.35 GeV,K=1.007) xl.( 
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FIG. 14: (Color online). Flux- integrated single differential 
cross section per target neutron for the CCQE process. 
The measured values are shown as points with the shape 
error as shaded bars. Calculations from the nuance RFG 
model with different assumptions for the model parameters 
are shown as histograms. Numerical values are provided in 
Table HXl in the Appendix. 



reeled for detector resolution effects only ( Section lIVD I) . 
Thus, this result is the most model-independent mea- 
surement of this process possible with the MiniBooNE 
detector. No requirements on the nucleonic final state 
are used to define this process. The neutrino flux is an 
absolute prediction [19] and has not been adjusted based 
on measured processes in the MiniBooNE detector. 



V. RESULTS AND DISCUSSION 

A. CCQE flux-integrated double differential cross 
section 

The flux-integrated, double differential cross section 
per neutron, dcos9 ' ^t^ CCQE process is ex- 

tracted as described in Section IIVDI and is shown in 
Figure [13] for the kinematic range, — 1 < cos 6*^ < +1, 
0.2 < r^(GeV) < 2.0. The errors, for outside of this 
range, are too large to allow a measurement. Also, bins 
with low event population near or outside of the kine- 
matic edge of the distribution (corresponding to large 
El,) do not allow for a measurement and are shown as 
zero in the plot. The numerical values for this double 
differential cross section are provided in Table IVII in the 
Appendix. 

The flux-integrated CCQE total cross section, ob- 
tained by integrating the double differential cross section 
(over -1 < cos^p < -1-1, < r^(GeV) < c»), is mea- 
sured to be 9.429 x lO^'^^ cm^. The total normalization 
error on this measurement is 10.7%. 

The kinematic quantities, and cos 9^, have been cor- 



B. Flux-integrated single differential cross section 

The flux-integrated, single differential cross section per 
neutron, -jM—, has also been measured and is shown 

in Figure. [Ml The quantity Qq^ is deflned in Eq. [2| 
and depends only on the (unfolded) quantities and 
cos9f^. It should be noted that the efficiency for events 
with < 200 MeV is not zero because of difference 
between reconstructed and unfolded T^. The calculation 
of efficiency for these (low-Qg^) events depends only on 
the model of the detector response, not on an interaction 
model and the associated uncertainty is propagated to 
the reported results. 

In addition to the experimental result. Figure [TD also 
shows the prediction for the CCQE process from the nu- 
ance simulation with three different sets of parameters 
in the underlying RFG model. The predictions are ab- 
solutely normalized and have been integrated over the 
MiniBooNE flux. The RFG model is plotted assum- 
ing both the world-avejaged CCQE parameters {Ma = 
1.03 GeV, K = 1.000) ;9fl and the CCQE parameters ex- 
tracted from this analysis (M^ = 1.35 GeV, k = 1.007) 
in a shape-only flt. The model using the world-averaged 
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FIG. 15: (Color online). Flux-unfolded MiniBooNE z/^, CCQE 
cross section per neutron as a function of neutrino energy. In 
(a), shape errors are shown as shaded boxes along with the 
total errors as bars. In (b), a larger energy range is shown 
along with results from the LSND 56] and NOMAD ;10j ex- 
periments. Also shown are predictions from the nuance sim- 
ulation for an RFG model with two different parameter vari- 
ations and for scattering from free nucleons with the world- 
average Ma value. Numerical values are provided in Table 1X1 
in the Appendix. 



CCQE parameters underpredicts the measured differen- 
tial cross section values by 20 — 30%, while the model 
using the CCQE parameters extracted from this shape 
analysis are within k, 8% of the data, consistent within 
the normalization error (w 10%). To further illustrate 
this, the model calculation with the CCQE parameters 
from this analysis scaled by 1.08 is also plotted and shown 
to be in good agreement with the data. 



C. Flux-unfolded CCQE cross section as a function 
of neutrino energy 

The flux-unfolded CCQE cross section per neutron, 
a[E'^^'^^^\, as a function of the true neutrino energy, 
j^QE,RFG ^ is shown in Figure [15] These numerical values 
are tabulated in Table |X] in the Appendix. The quantity 
j^QE,RFG jg (model-dependent) estimate of the neu- 
trino energy obtained after correcting for both detector 
and nuclear model resolution effects. These results de- 
pend on the details of the nuclear model used for the cal- 
culation. The dependence is only weak in the peak of the 
flux distribution but becomes strong for Ei, < 0.5 GeV 
and E„ > 1.2 GeV, i.e., in the "tails" of the flux distri- 
bution. 

In Figure [121 the data are compared with the nuance 
implementation of the RFG model with the world average 
parameter values, {M'f = 1.03 GeV, k = 1.000) and 
with the parameters extracted from this work (M^* = 
1.35 GeV, K = 1.007). These are absolute predictions 
from the model (not scaled or renormalized). At the 



source 


normalization error (%) 


neutrino flux prediction 


8.66 


background cross sections 


4.32 


detector model 


4.60 


kinematic unfolding procedure 


0.60 


statistics 


0.26 


total 


10.7 



TABLE IV: Contribution to the total normalization uncer- 
tainty from each of the various systematic error categories. 



average energy of the MiniBooNE flux (w 800 MeV) , the 
extracted cross section is 30% larger than the RFG 
model prediction with world average parameter values. 
The RFG model, with parameter values extracted from 
the shape- only fit to this data better reproduces the data 
over the entire measured energy range. 

Figure [TSTb) shows these CCQE results together with 
those from the LSND |56i and NOMAD [T^ experiments. 
It is interesting to note that the NOMAD results are bet- 
ter described with the world-average AI"^ and k values. 
Also shown for comparison in Fig. I15f b) is the predicted 
cross section assuming the CCQE interaction occurs on 
free nucleons with the world-average Ma value. The cross 
sections reported here exceed the free nucleon value for 
E,, above 0.7 GeV. 



D. Error Summary 



As described in Section HV El (correlated) systematic 
and statistical errors are propagated to the final results. 
These errors are separated into normalization and shape 
uncertainties. The contributions from each error source 
on the total normalization uncertainty are summarized 
in Table IIVI As is evident, the neutrino flux uncer- 
tainty dominates the overall normalization error on the 
extracted CCQE cross sections. However, the uncer- 
tainty on the flux prediction is a smaller contribution 
to the shape error on the cross sections. This can be 
seen in Figure [TBI which shows the contribution from the 
four major sources to the shape error on the total (flux- 
unfolded) cross section. 

The detector model uncertainty dominates the shape 
error, especially at low and high energies. This is because 
errors in the detector response (mainly via uncertain- 
ties in visible photon processes) will result in errors on 
the reconstructed energy. These errors grow in the tails 
of the neutrino flux distribution due to feed-down from 
events in the fiux peak. This type of measurement usu- 
ally has large errors due to non-negligible uncertainties 
in the CCl7r+ background predictions. In this measure- 
ment, that error is reduced through direct measurement 
of the CCl7r+ background. However, this error is not 
completely eliminated due to the residual uncertainty on 
the rate of intranuclear pion absorption that is included. 
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FIG. 16: Fractional shape error on the MiniBooNE Ufj, CCQE 
flux-unfolded cross section separated into major components. 
The overall normalization error of 10.7% is not shown. 

This uncertainty is not as important for the measurement 
of the CCQE cross section measurement as a function of 
energy but is a large contribution to the error at low Qq^ 
in the differential distributions. 

The unfolding error is small in the region of the flux 
peak but grows in the high- and low-energy region be- 
cause of the uncertainty in the feed-down from other 
energy bins, similar to that described for the detector 
model. 



VI. CONCLUSIONS 

In this work, we report measurements of absolute cross 
sections for the CCQE interaction using high-statistics 
samples of interactions on carbon. These include the 
first measurement of the double differential cross section, 

.2 

-ppA-^s-, measurement of the single differential cross 

dTf^d cos 6^ ' ^ 

section, .i?? , and the flux- unfolded CCQE cross sec- 

tion, ct[S^^'-«^'^]. The double differential cross section 
contains the most complete and model-independent infor- 
mation that is available from MiniBooNE on the CCQE 
process. It is the main result from this work and should 
be used as the preferred choice for comparison to theo- 
retical models of CCQE interactions on nuclear targets. 

The reported cross section is signiflcantly larger (« 
30% at the flux average energy) than what is commonly 
assumed for this process assuming a relativistic Fermi 
Gas model (RFC) and the world- average value for the ax- 
ial mass, Ma = 1.03 GeV In addition, the Qq^ dis- 
tribution of this data shows a signiflcant excess of events 
over this expectation at higher Qq^ even if the data is 
normalized to the prediction over all Qqe- This leads 
to an extracted axial mass from a "shape-only" flt of the 
Q^^ distribution of M^f = 1.35±0.17 GeV, significantly 
higher than the historical world-average value. 

These two observations, unexpectedly large values for 



the extracted cross section and , are experimen- 
tally separate. However, within the model prediction, 
a larger value for Ma implies a larger cross section be- 
cause the CCQE cross section increases approximately 
linearly with AIa- The predicted CCQE cross section 
with this higher value of Ma agrees with the measure- 
ment within the normalization error of the experiment 
(~ 10%). While this may be simply a coincidence, it is 
important to note. 

In recent years, there has been significant effort to im- 
prove the theoretical description of the CCQE interaction 
on nuclear targets [H, [13; however, there seems to be 
no simple explanation for both the higher cross section 
and the harder Qq^ distribution of events (resulting in a 
larger Mf) as evidenced by the MiniBooNE CCQE data. 
Nuclear effects can have some impact on the measured 
M'^, but it is not obvious that they are large enough. 
Also, it is expected that such effects should reduce the 
cross section, not increase it. This can be be seen in 
Fig. [T51 where the cross section for the CCQE interaction 
on free nucleons is compared to that from bound nu- 
cleons (in an RFC model). Note the reduction in cross 
section from free to bound nucleons. It is interesting that 
the MiniBooNE measurement is also larger than this free 
nucleon value (at least at higher energies). This may in- 
dicate a signiflcant contribution from neglected mecha- 
nisms for CCQE-like scattering from a nucleus such as 
multi- nucleon processes (for example, Ref. [l3l)- This 
may explain both the higher cross section and the harder 
spectrum, but has not yet been explicitly tested. It 
may also be relevant for the difference between these re- 
sults and those of NOMAD (or other experiments) where 
the observation of recoil nucleons enter the deflnition of 
a CCQE event. An important test for such models will 
be their ability to accurately reproduce the MiniBooNE 
double differential cross sections at least as well as the 
RFG model assuming a higher axial mass value. 

As yet, there is no easily recognized solution to explain 
the difference between the CCQE cross sections mea- 
sured in MiniBooNE at lower neutrino energy (i?^ < 2 
GeV) and the NOMAD results at higher neutrino ener- 
gies {El, > 3 GeV). Model- independent measurements of 
the CCQE cross section anticipated from SciBooNE [57| - 
MicroBooNE and MINERvA ^ as well as the 
T2K @ and NOvA ^ near detectors running with 
2 < < 20 GeV, will be important to help resolve 
these results. 
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Appendix A: Tabulation of results 

This appendix contains tables of numerical values cor- 
responding to various plots appearing in the main body 
of the paper. In addition, the effective cross section for 
the CCQE-like background to the CCQE measurement 
is tabulated. These tables are also available via the Mini- 
BooNE website 

1. Predicted Uf^ flux 

Table |V] lists the predicted flux (Figure [2]) at 
the MiniBooNE detector in 50 MeV-wide neutrino en- 
ergy bins. The flux is normalized to protons on target 
(POT). The mean energy is 788 MeV and the integrated 
flux over the energy range (0.0 < < 3.0 GeV) is 
5.16 X 10-1° i/^/POT/cm^. For this analysis, the to- 
tal POT collected is 5.58 x 10^° yielding an integrated 
flux of 2.88 X 10" ly^/cm^. 

2. CCQE flux-integrated double differential cross 
section 

Table I VII contains the flux-integrated CCQE dou- 
ble differential cross section values ( .^^ ^ „ ) in bins of 

^ dTf^ a cos 

muon energy, T^, and cosine of the muon scattering an- 
gle with respect to the incoming neutrino direction (in 
the lab frame), cos0^. These values correspond to the 
plot of Figure [TS] The integrated value over the region ( 
-1 < cos 6*^ < -1-1 and < < oo) is 9.429x10^39 ^,^^2 
The total normalization uncertainty is 10.7%. Table IVlTl 
present an analagous summary of the shape error for each 
bin. 



3. CCQE-like backgrounds 

As explained in Sect. lIVDl the CCl7r+ interaction with 
intranuclear pion absorption forms a "CCQE-like" back- 



ground in that the final state is indistinguishable from 
the CCQE signal in MiniBooNE. These events originate 
from the CCItt"*" interaction but contain 1 muon and 
no pions in the final state. In the main analysis, this 
background is subtracted to obtain the CCQE observ- 
ables. In order to facilitate comparisons with models 
(or other analyses) that consider all CCQE-like events as 
CCQE signal, the effective double differential cross sec- 
tion for the CCItt^ interaction with intranuclear pion 
absorption is presented in Table I Villi These values are 
determined from the NUANCE-event generator corrected 
to reproduce the MiniBooNE 3-subevent sample and are 
calculated using Eq. [3] with {dj — bj ) replaced by b'j , the 
number of CCQE-like background events. A CCQE-like 
cross section may be obtained by adding these numbers 
(Table IVml with those from Table |Vll 
4. CCQE flux-integrated single differential cross 
section 



Table lIXl contains the flux-integrated CCQE single dif- 
ferential cross section ( ) in bins of Qqe- Qqe 

defined in Eq. [2j The shape error and CCQE-like back- 
ground prediction is also reported. The corresponding 
plot is shown in Figure 1141 



5. Flux unfolded CCQE cross section 



TablelXlcontains the flux-unfolded CCQE cross sec- 
tion values bins of EQ^'^^^. eQ^'^^^^ 
is as deflned in Eq. [T] The shape and total errors along 
with the CCQE-like background are also reported. The 
corresponding plot is shown in Figure 1151 
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flux 
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Iy'f_i flux 


El, bin 


flux 


(GeV) 


(i/^/POT/GeV/cm^) 


(GeV) 


(i/^/POT/GeV/cm^) 


(GeV) 


(i^^/POT/GeV/cm^) 


0.00-0.05 


4.54 


X 


10-" 


1.00-1.05 


3.35 


X 


10-10 


2.00-2.05 
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X 
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1.71 
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1.05-1.10 
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10-10 
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0.25-0.30 
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2.14 
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2.25-2.30 
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1.30-1.35 
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1.70-1.75 
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X 
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2.28 
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3.20 


X 
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TABLE V: Predicted Uf, flux at the MiniBooNE detector. 



0^T^{GcV) 0.2,0.3 0.3,0.4 0.4,0.5 0.5,0.6 0.6,0.7 0.7,0.8 0.8,0.9 0.9,1.0 1.0,1.1 1.1,1.2 1.2,1.3 1.3,1.4 1.4,1.5 1.5,1.6 1.6,1.7 1.7,1.8 1.8,1.9 1.9,2.0 



+0.9,-1-1.0 


190.0 


326.5 


539.2 


901.8 


1288 


1633 


1857 


1874 


1803 


1636 


1354 


1047 


794.0 


687.9 


494.3 


372.5 278.3 227. 


+0.8, +0.9 


401.9 


780.6 


1258 


1714 


2084 


2100 


2035 


1620 


1118 


783.6 


451.9 


239.4 


116.4 


73.07 


41.67 


36.55 — — 


+0.7, +0.8 


553.6 


981.1 


1501 


1884 


1847 


1629 


1203 


723.8 


359.8 


156.2 


66.90 


26.87 


1.527 


19.50 






+0.6, +0.7 


681.9 


1222 


1546 


1738 


1365 


909.6 


526.7 


222.8 


81.65 


35.61 


11.36 


0.131 










+0.5, +0.6 


765.6 


1233 


1495 


1289 


872.2 


392.3 


157.5 


49.23 


9.241 


1.229 


4.162 












+0.4, +0.5 


871.9 


1279 


1301 


989.9 


469.1 


147.4 


45.02 


12.44 


1.012 
















+0.3, +0.4 


910.2 


1157 


1054 


628.8 


231.0 


57.95 


10.69 




















+0.2, +0.3 


992.3 


1148 


850.0 


394.4 


105.0 


16.96 


10.93 




















+0.1, +0.2 


1007 


970.2 


547.9 


201.5 


36.51 


0.844 






















0.0, +0.1 


1003 


813.1 


404.9 


92.93 


11.63 
























-0.1, 0.0 


919.3 


686.6 


272.3 


40.63 


2.176 
























-0.2,-0.1 


891.8 


503.3 


134.7 


10.92 


0.071 
























-0.3,-0.2 


857.5 


401.6 


79.10 


1.947 


























-0.4,-0.3 


778.1 


292.1 


33.69 




























-0.5,-0.4 


692.3 


202.2 


17.42 




























-0.6,-0.5 


600.2 


135.2 


3.624 




























-0.7,-0.6 


497.6 


85.80 


0.164 




























-0.8,-0.7 


418.3 


44.84 






























-0.9,-0.8 


348.7 


25.82 






























-1.0,-0.9 


289.2 


15.18 































TABLE VI: The MiniBooNE CCQE flux-integrated double differential cross section in units of 10-*i cm2/GeV in 0.1 GeV 
bins of (columns) and 0.1 bins of cos 6^ (rows). 



1189, 133 (2009); M. E. Dorman, FERMILAB-THESIS- 
2008-72. 

[16] J. E. Amaro et at, Phys. Rev. C 71, 015501 (2005); 
Phys. Rev. C 75, 034613 (2007); A. M. Ankowski et 
al, Phys. Rev. C 77, 044311 (2008); O. Benhar et 



al., Phys. Rev. D 72, 053005 (2005); arXiv:0903.2329J 
[hep-ph]; A. Butkevich et al, Phys. Rev. C 76, 045502 
(2007); Phys. Rev. C 80, 014610 (2009); N. Jachowicz 
et al, Phys. Rev. C 73, 024607 (2006); T. Leitner et 
al, Phys. Rev. C 73, 065502 (2006); Phys. Rev. C 79, 
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0.2,0.3 0.3,0.4 0.4,0.5 0.5,0.6 0.6,0.7 0.7,0.8 0.8,0.9 0.9,1.0 1.0,1.1 1.1,1.2 1.2,1.3 1.3,1.4 1.4,1.5 1.5,1.6 1.6,1.7 1.7,1.8 1.8,1.9 1.9,2.0 

684.3 1071 1378 1664 1883 2193 2558 3037 3390 3320 3037 3110 2942 2424 2586 2653 3254 3838 
905.0 1352 1754 2009 2222 2334 2711 2870 2454 1880 1391 1036 758.7 544.3 505.5 359.6 — — 
1134 1557 1781 1845 1769 1823 1873 1464 963.8 601.6 339.6 184.1 170.1 230.6 — — — — 

1435 1455 1581 1648 1791 1513 1068 598.2 267.2 155.1 69.28 89.01 — — — — — — 

1380 1372 1434 1370 1201 870.2 432.3 162.2 71.88 49.10 54.01 — — — — — — — 

1477 1273 1365 1369 1021 475.5 161.6 55.58 16.32 — — — — — — — — — 

1267 1154 1155 965.3 574.7 149.2 53.26 — — — — — — — — — — — 

1293 1105 1041 742.5 250.6 77.66 110.3 — — — — — — — — — — — 

1351 1246 1048 415.1 114.3 41.02 — — — — — — — — — — — — 

1090 1078 695.5 238.2 45.96 — — — — — — — — — — — — — 

980.4 783.6 515.7 114.6 20.92 — — — — — — — — — — — — — 

917.7 746.9 337.5 50.92 3.422 — — — — — — — — — — — — — 

922.7 586.4 215.6 55.88 — — — — — — — — — — — — — — 

698.0 553.3 135.3 — — — — — — — — — — — — — — — 

596.9 482.6 57.73 — — — — — — — — — — — — — — — 

520.8 360.7 34.63 — — — — — — — — — — — — — — — 

450.2 236.6 31.22 — — — — — — — — — — — — — — — 

408.8 184.4 — — — — — — — — — — — — — — — — 

339.7 107.6 — — — — — — — — — — — — — — — — 

349.8 63.32 — — — — — — — — — — — — — — — — 



TABLE VII: Shape uncertainty on the MiniBooNE i/^ CCQE flux-integrated double differential cross section in units of 
10^*^ cm^/GeV corresponding to Table IVll The total normalization error is 10.7%. 



g Tj,(GgV) 0.2,0.3 0.3,0.4 0.4,0.5 0.5,0.6 0.6,0.7 0.7,0.8 0.8,0.9 0.9,1.0 1.0,1.1 1.1,1.2 1.2,1.3 1.3,1.4 1.4,1.5 1.5,1.6 1.6,1.7 1.7,1.8 1.8,1.9 1.9,2.0 



-1-0.9,-1-1.0 


83.6 


199.8 


285.3 


364.2 


391.1 


403.7 


384.3 


349.2 


301.4 


232.7 


179.2 


136.1 


102.0 


90.73 


76.55 


52.36 41.47 54.50 


+0.8,-1-0.9 


111.6 


257.4 


351.0 


364.3 


353.2 


288.9 


233.8 


169.5 


106.6 


59.81 


31.21 


20.89 


10.10 


6.008 


2.376 


2.859 — — 


+0.7,-1-0.8 


118.4 


270.4 


312.6 


280.3 


211.7 


135.7 


81.47 


40.97 


21.56 


9.247 


3.284 


0.875 


0.057 








+0.6, +0.7 


118.9 


260.0 


252.8 


183.4 


101.8 


52.52 


19.75 


7.978 


2.716 


0.281 














+0.5, +0.6 


109.0 


215.2 


181.4 


104.6 


41.87 


16.33 


3.643 


0.492 


0.004 
















+0.4, +0.5 


109.2 


182.0 


122.4 


51.26 


19.76 


4.193 


0.183 




















+0.3, +0.4 


104.0 


140.2 


73.71 


24.54 


4.613 


0.151 


0.002 




















+0.2, +0.3 


93.84 


107.6 


48.56 


10.78 


0.812 
























+0.1, +0.2 


76.55 


80.94 


29.02 


3.049 


0.030 
























0.0, +0.1 


67.81 


52.89 


13.71 


0.392 


























-0.1, 0.0 


58.91 


37.46 


5.565 


0.011 


























-0.2,-0.1 


50.47 


22.49 


1.048 




























-0.3,-0.2 


39.03 


12.58 


0.118 




























-0.4,-0.3 


32.41 


7.575 


0.061 




























-0.5,-0.4 


25.72 


2.529 


0.080 




























-0.6,-0.5 


16.78 


1.063 


0.009 




























-0.7,-0.6 


9.963 


0.280 


0.002 




























-0.8,-0.7 


5.005 


0.244 






























-0.9,-0.8 


4.877 


0.067 






























-1.0,-0.9 


3.092 


0.013 































TABLE VIII: The predicted CCQE-like background flux-integrated double differential cross section in units of 
10""^ cmVCeV corresponding to Table IVll 
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Q%E (GeV^ 



shape error (cmVCeV^) CCQE-like bkgd (cmVCeV^ 



0.00-0.05 


7.681 


X 


1 n— 39 

10 


1 


.493 


X 


1 n — 39 
10 


3.876 


X 


1 n — 39 

10 


0.05-0.10 


1.457 


X 


10-3* 


1 


.180 


X 


^ r>— 39 
10 


3.961 


X 


1 n — 39 
10 


0.10-0.15 


1.684 


X 


10-3* 


9, 


.720 


X 


T ri — 40 
10 


3.671 


X 


^ r> — 39 

10 


0.15-0.20 


1.703 


X 


10-3* 


8, 


.216 


X 


1 n — 40 
10 


3.064 


X 


T ri— 39 
10 


0.20-0.25 


1.589 


X 


1 n— 38 
10 


5 


.134 


X 


T ri — 40 
10 


2.522 
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1 n — 39 
10 


0.25-0.30 


1.449 


X 


10-3* 
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.983 
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T ri — 40 
10 


2.040 


X 


^ r> — 39 
10 


0.30-0.35 


1.329 
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10-3* 


3 


.386 


X 


T ri — 40 
10 


1.633 


X 


1 n — 39 
10 


0.35-0.40 


1.172 


X 


10-3* 


2, 


.629 


X 


T ri — 40 
10 


1.290 


X 


1 n — 39 
10 


0.40-0.45 


1.030 


X 


10-3* 


2, 


.457 


X 


T ri — 40 
10 


1.018 


X 


1 r> — 39 
10 
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1 n — 39 

10 
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.975 
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T ri— 40 

10 


7.874 
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T ri — 40 
10 
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T ri — 40 
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.212 


X 


10-''° 
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.442 
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.504 


X 
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X 
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TABLE IX: The MiniBooNE CCQE flux-integrated single differential cross section, errors, and predicted CCQE-like 
background in bins of Qqe- The total normalization error is 10.7%. 



(GeV) (7 (cm^) shape error (cm^) total error (cm^) CCQE-like bkgd (cm^ 



0.40-0.45 


7, 
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10" 


-39 


1.997 


X 


10" 


39 


1 


.997 


X 
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-39 


1, 


.731 


X 
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-39 
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.261 


X 


10" 


-39 


1.455 


X 
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1 


.532 


X 
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-39 


1, 


.865 


X 


10" 


-39 


0.50-0.55 


8, 


.809 


X 


10- 


-39 


1.169 


X 


10" 


-39 


1 


.330 


X 
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-39 


1, 


.951 


X 


10" 


-39 


0.55-0.60 


9, 


.530 


X 


10" 


-39 


9.537 


X 
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-40 


1 


.209 


X 
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-39 


1, 


.978 


X 


10" 
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1, 


.013 


X 


10" 


-38 


7.575 


X 
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-40 


1 


.124 


X 


10" 


-39 


1, 


.941 


X 


10" 


-39 


0.65-0.70 


1, 


.071 


X 


10" 


-38 


6.000 


X 


10" 


-40 


1 


.089 


X 


10" 


-39 


1, 


.878 


X 


10" 


-39 
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1, 


.111 


X 


10" 


-38 


4.496 


X 


10" 


-40 


1 


.065 


X 


10" 


-39 


1, 


.770 


X 


10" 


-39 
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1, 


.155 


X 
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-38 
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X 
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.078 
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.672 
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.278 


X 


10" 
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X 
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X 
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7, 
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X 
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-40 
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.236 


X 


10" 


-38 
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X 


10" 


-39 


2 


.613 


X 


10" 


-39 


4, 


.822 


X 


10" 
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TABLE X: The MiniBooNE CCQE flux-unfolded cross section, errors, and predicted CCQE-like background in bins of 

j^QE,RFG 
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